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Observational evidence and theoretical motivation for dark matter are presented 
and connections to the CMB and BBN are made. Problems for baryonic and 
neutrino dark matter are summarized. Emphasis is placed on the prospects for 
supersymmetric dark matter. 



1. Lecture 1 

The nature and identity of the dark matter of the Universe is one of the 
most challenging problems facing modern cosmology. The problem is a 
long-standing one, going back to early observations of mass-to-light ratios 
by ZwickjW Given the distribution (by number) of galaxies with total lu- 
minosity L, 4>{L), one can compute the mean luminosity density of galaxies 

C = J L<p(L)dL (1) 

which is determined to bJ^ 

£ ~ 2 ± 0.2 x l0 s h o L Q Mpc- 3 (2) 

where Lq = 3.8 x 10 33 erg s" 1 is the solar luminosity. In the absence 
of a cosmological constant, one can define a critical energy density, p c — 
3H 2 /8ttGn = 1-88 x l0~ 29 h o 2 g cm~ 3 , such that p = p c for three-space 
curvature k — 0, where the present value of the Hubble parameter has 
been defined by H a = 100/i o km Mpc -1 s _1 . We can now define a critical 
mass-to-light ratio is given by 

(M/L) c = Pc /C ~ 1390MA/O/LQ) (3) 



'Summary of lectures given at the Theoretical Advanced Study Institute in Elementary 

Particle Physics at the University of Colorado at Boulder - June 2-28, 2002. 
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which can be used to determine the cosmological density parameter 



Mass-to-light ratios are, however, strongly dependent on the distance 
scale on which they are determined^. In the solar neighborhood M/L ~ 
2 ± 1 (in solar units), yielding values of Q m of only ~ .001. In the bright 
central parts of galaxies, M/L ~ (10 — 20)h o so that tt m ~ 0.01. On larger 
scales, that of binaries and small groups of galaxies, M/L ~ (60 — 180) h a 
and D, m ~ 0.1. On even larger scales, that of clusters, M/L may be as 
large as (200 — 500)/i o giving fi m ~ 0.3. This progression in M/L seems 
to have halted, as even on the largest scales observed today, mass-to-light 
ratios imply values of fl m < 0.3 — 0.4. Thus when one considers the scale 
of galaxies (and their halos) and larger, the presence of dark matter (and 
as we shall see, non-baryonic dark matter) is required. 

Direct observational evidence for dark matter is found from a variety 
of sources. On the scale of galactic halos, the observed flatness of the 
rotation curves of spiral galaxies is a clear indicator for dark matter. There 
is also evidence for dark matter in elliptical galaxies, as well as clusters of 
galaxies coming from the X-ray observations of these objects. Also, direct 
evidence has been obtained through the study of gravitational lenses. On 
the theoretical side, we predict the presence of dark matter (or dark energy) 
because 1) it is a strong prediction of most inflation models (and there is at 
present no good alternative to inflation) and 2) our current understanding 
of galaxy formation requires substantial amounts of dark matter to account 
for the growth of density fluctuations. One can also make a strong case for 
the existence of non-baryonic dark matter in particular. The recurrent 
problem with baryonic dark matter is that not only is it very difficult to 
hide baryons, but given the amount of dark matter required on large scales, 
there is a direct conflict with primordial nucleosynthesis if all of the dark 
matter is baryonic. In this first lecture, I will review the observational and 
theoretical evidence supporting the existence of dark matter. 

1.1. Observational Evidence 

Assuming that galaxies are in virial equilibrium, one expects that one can 
relate the mass at a given distance r, from the center of a galaxy to its 
rotational velocity by 



fi m = JL = (M/L)/{M/L) C 



Pc 



(4) 



M(r) oc v 2 r/G N 



(5) 
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Figure 1. Synthetic rotation curvJSfor galaxies with (M) = —21.2. The dotted curve 
shows the disk contribution, whereas the dashed curve shows the halo contribution. 



The rotational velocity, v, is measured^ - by observing 21 cm emission lines 
in HI regions (neutral hydrogen) beyond the point where most of the light 
in the galaxy ceases. A subset of a compilation^ of nearly 1000 rotation 
curves of spiral galaxies is shown in Fig. ^ The subset shown is restricted 
to a narrow range in brightness, but is characteristic for a wide range of 
spiral galaxies. Shown is the rotational velocity as a function of r in units 
of the optical radius. If the bulk of the mass is associated with light, then 
beyond the point where most of the light stops, M would be constant and 
v 2 oc 1/r. This is not the case, as the rotation curves appear to be flat, 
i.e., v ~ constant outside the core of the galaxy. This implies that M oc r 
beyond the point where the light stops. This is one of the strongest pieces of 
evidence for the existence of dark matter. Velocity measurements indicate 
dark matter in elliptical galaxies as welP. 

Galactic rotation curves are not the only observational indication for the 
existence of dark matter. X-ray emitting hot gas in elliptical galaxies also 
provides an important piece of evidence for dark matter. A particularly 
striking example is that of the large elliptical M87. Detailed profiles of the 
temperature and density of the hot X-ray emitting gas have been mapped 
oulrJ Assuming hydrostatic equilibrium, these measurements allow one 
to determine the overall mass distribution in the galaxy necessary to bind 
the hot gas. Based on an isothermal model with temperature kT = 3keV 
(which leads to a conservative estimate of the total mass), Fabricant and 
GorensteirPI predicted that the total mass out to a radial distance of 392 
kpc is 5.7 x 10 13 Af Q , whereas the mass in the hot gas is only 2.8 x 1O 12 M 
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or only 5% of the total. The visible mass is expected to contribute only 1% 
of the total. 

M87 is not the only example of an elliptical galaxy in which X-ray emit- 
ting hot gas is observed to indicate the presence of dark matter. X-ray ob- 
servations have shown that the total mass associated with elliptical galaxies 
is considerably larger than the luminous component in many examples of 
varying morphological typeJ3. Mass-to-light ratios for these systems vary, 
with most being larger than 30/io and some ranging as high as ~ 200/io- 

In addition, similar inferences pertaining to the existence of dark matter 
can be made from the X-ray emission from small groups of galaxiea^l. On 
these scales, mass-to-light ratios are typically > lOO/io and detailed studies 
have shown that the ba ryo n faction in these systems is rather small. Fur- 
thermore, it was argued^ that cluster baryon fractions should not differ 
from the Universal value given by Sls/Slm- Using an estimate of tts = 0.04 
from big bang nucleosynthesis (BBN, see below), and baryon fractions rang- 
ing from 0.1 to 0.3, one would obtain an estimate for the total matter density 
of 0.13 - 0.4. 

Another piece of evidence on large scales is available from gravitational 
lensin^^. The systematic lensing of the roughly 150,000 galaxies per deg 2 
at redshifts between z = 1 — 3 into arcs and arclets allow one to trace 
the matter distribution in a foreground cluster. Lensing observations can 
be categorized as either strong (multiple images) or weak (single images). 
Both require the presence of a dominant dark matter component. 

Strong lensing is particularly adept in testing the overall geometry of 
the Univers e ! I I While a cluster which provides multiple lenses of a single 
background galaxy (at known redshift) is useful for determining the total 
cluster mass, when several background galaxies are lensed, it is possible to 
constrain the values of fi m and SI a "I The recent results of 021 show a 
degeneracy in the Q m - 51a plane. Nevertheless, the allowed region is offset 
from similar types of degeneracies found in supernovae searches and the 
CMB (see below). Indeed, these lensing results are much more constraining 
in f2 m than the other techniques, though a residual uncertainty of about 
30 % persists. While in principle, these results find that any value (from 
to 1) is possible for Q m < 0.5 for low values of Ha and O m < 0.4 for 
higher values of Ha (k, 0.6). 

Weak lensing of galaxies by galaxies can (on a statistical basis) also 
probe the nature of galactic halos. Recent studies based on weak lensing 
data in dica te that galactic halos may be far more extended than previously 
though (radii larger than 200 h 1 kpc). These results also imply a sub- 
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stantial contribution to Q m (of order 0.1-0.2) on this scale. On larger scales, 
using many cluster lenses enables one to estimate ft ~ 0.3 ^M. Another use 
of weak lensing statistics is t o determine the evolution of cosmic shear and 
hence an estimate of Vt m Finally, there exist a number of examples of 
dark clusters, ie., lenses with no observable counterparlP^I The contribu- 
tion of these objects (if they are robust) to fl m is not clear. For a recent 
review of weak lensing see 

Finally, on very large scales, it is possible to get an estimate of fl m 
from the distribution of peculiar velocities of galaxies and clusters. On 
scales, A, where perturbations, 6, are still small, peculiar velocities can be 
expressed^ as v ~ HXSQ 1 ^ 6 . On these scales, older measurements of the 
peculiar velocity field from the IRAS galaxy catalogue indicate that indeed 
f2 is clos e to unity" 1 . Some of the new data indicates a low erva lue in the 
range^2 0.2 - 0.5, but does not conclusively exclude Q m = 1 1=^! 

The above discussion of observational evidence for dark matter pertains 
largely to the overall matter density of the Universe f2 m . (F or a comprehen- 
sive review of determinations of the matter density, see ref. 1=21 ) However, 
the matter density and the overall curvature are not related one-to-one. The 
expansion rate of the Universe in the standard FRW model is expressed by 
the Friedmann equation 

2 &_ 8ttG nP _ k_ A 

R 2 3 R? 3 [ 1 

where R(t) is the cosmological scale factor, k is the three-space curvature 
constant (k = 0, +1,— 1 for a spatially flat, closed or open Universe), and 
A is the cosmological constant. The Friedmann equation can be rewritten 
as 

(n - m 2 = ^ (7) 

so that k = 0, +1, —1 corresponds to £1 = 1, f2 > 1 and f2 < 1. However, 
the value of Q appearing in Eq. (JTJ represents the sum Q = Sl m + of 
contributions from the matter density (Cl m ) and the cosmological constant 
(fi A = A/3H 2 ). 

There has been a great deal of progress in the last several years concern- 
ing the determination of both O m and Q\. Cosmic Microwave Background 
(CMB) anisotropy experiments have been able to determine the curvature 
(i.e. the sum of 17 m and Q\) to with in about 10%, while observations 
of type la supernovae at high redshift provide information on a (nearly) 
orthogonal combination of the two density parameters. 
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The CMB is of course deeply rooted in the development and verification 
of the big bang modeP^l. Indeed, it was the formulation of BBN that led 
to the prediction of the microwave background. The argument is rather 
simple. BBN requires temperatures greater than 100 keV, which according 
to the standard model time-temperature relation, t s T^ c y = 2A/y/~N, where 
N is the number of relativistic degrees of freedom at temperature T, and 
corresponds to timescales less than about 200 s. The typical cross section 
for the first link in the nucleosynthetic chain is 



for nucleosynthesis to begin. The density in baryons today is known ap- 
proximately from the density of visible matter to be Ub ~ 10 -7 cm -3 and 
since we know that that the density n scales as R~ 3 ~ T 3 , the temperature 
today must be 



thus linking two of the most important tests of the Big Bang theory. 

Microwave background anisotropy measurements have mad e t rem en- 
dous advances in the last few years. The power 

spectruirPHEZIlD 

has 

been measured relatively accurately out to multipole moments corre- 
sponding to £ ~ 1000. The details of this spectrum enable one to 
make accurate predicti ons of a large number of fundamental cosmologi- 
cal parameter J^7 | 29 | 30 | 31 | 32 (55| ^ n exanl pi e f these results as found by a 
recent frequentist analysiJSl is shown in Fig. [5] 

Of particular interest to us here is the CMB determination of the total 
density, Sl to t, as well as the matter density O m . The results of recent CMB 
anisotropy measurements are summarized in Table 1. As one can see, there 
is strong evidence that the Universe is flat or very close to it. Furthermore, 
the matter density is very consistent with the observational determinations 
discussed above and the baryon density, as we will see below, is consistent 
with the BBN production of D/H and its abundance in quasar absorption 
systems. 

The discrepancy between the CMB value of fi m and VIb is sign that 
non-baryonic matter (dark matter) is required. Furthermore, the apparent 
discrepancy between the CMB value of f2 to t an d ^rm though not conclusive 
on its own, is a sign that a contribution from the vacuum energy density or 




(8) 



(9) 



crvt 



T = (n Bo /n) 1 / 3 T B BN ~ 10K 



(10) 
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Figure 2. Ax 2 calculated with the MAXIMA-1 and COBE data as a function of pa- 
rameter value. Solid blue circles show grid points in parameter space, and the green 
lines were obtained by interpolating between grid points. The parameter values where 
the green line intercepts th e re d dashed (dotted) line corresponds to the 68% (95%) 
frequentist confidence regional 



Table 1. Results from recent CMB anisotropy measurements. 





S^tot 




n B h 2 


BOOMERandl3 
MAXImJUI 
MAXIMA (freq.JSl 
DASP1 
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1.03 ±0.06 

0.891°;^ 

1.04 ±0.06 
0.99 ±0.12 
1.03 ±0.012 

1 lfi+ ' 24 

■ L - lo -0.30 


0.12 ±0.05 

n 1 7 +0.08 

25+ - 07 

0.14 ±0.04 

17 +0.08 
u ' ±( -0.06 
n 10 + O.O8 
u - iJ -0.05 


021 + 004 
0.0325 ±0.0125 
026+ 010 
022+ - 004 
022+ 015 
0.029 ±0.009 
u.uiy_ 007 



cosmological constant, is also required. The preferred region in the f2 m — SI a 
plane as determined by a frequentist analysis of MAXIMA data is shown 
in Fig. 
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Figure 3. Two-dimensional frequentist confidence regions in the (Hmi^a) plane. The 
red, orange and yellow regions correspond to the 68%, 95%, and 99% confidence regions 
respectively. The dashed black line corresponds to a flat universe, Q = f2 m + Oa = !• 



The presence or absence of a cosmological constant is a long standing 
problem in cosmology. To theorists, it is particularly offensive due to the 
necessary smallness of the constant. We know that the cosmological term is 
at most a factor of a few times larger than the current mass density. Thus 
from Eq. ©, we see that the dimensionless combination, GjyA < 10~ 121 . 
Nevertheless, even a small non-zero value for A could greatly affect the 
future history of the Universe: allowing open Universes to recollapse (if 
A < 0), or closed Universes to expand forever (if A > and sufficiently 
large). 

Another exciting development has been the use of type la supernovae, 
which now allow measurement of relative distances with 5% precision. In 
combination with Cepheid data from the HST key project on the distance 
scale, SNe results are the dominant contributor to the best modern value 
for H : 72kms _1 Mpc~ 1 ± 10%«*l Better still, the analysis of high-z SNe 
has allowed the first meaningful test of cosmological geometry to be carried 
out, as shown in Fig. 

These results can be contrasted with those from the CMB anisotropy 
measurements as in Fig. We are led to a seemingly conclusive picture. 
The Universe is nearly flat with Qtot — 1- However the density in matter 
makes up only 20-50% of this total, with the remainder in a cosmological 
constant or some other form of dark energy. 
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Figure 4. The type la supernova Hubble diagrarJSSl taken from The first panel 
shows that, for z < 1, the large-scale Hubble flow is indeed linear and uniform; the 
second panel shows an expanded scale, with the linear trend divided out, and with the 
redshift range extended to show how the Hubble law becomes nonlinear. Comparison 
with the prediction of Fricdmann-Lemaitre models appears to favor a vacuum-dominated 
universe. 

1.2. Theory 

Theoretically, there is no lack of support for the dark matter hypothesis. 
The standard big bang model including inflation almost requires fitot = 
1 1281. This can be seen from the following simple solution to the curvature 
problem. The simple and unfortunate fact that at present we do not even 
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Figure 5. Likelihood-based confidence contours^J over the plane f1yv (i- e - ^ 

u assuming 

w = — 1 vs Q m - The SNe la results very nearly constrain Cl v — fl m , whereas the results of 
CMB anisotropics (from the Boomerang 98 data) favor a flat model with Q v + Q m ~ 1. 
The intersection of these constraints is the most direct (but far from the only) piece of 
evidence favoring a flat model with f2 m ~ 0.3. 



know whether £1 is larger or smaller than one, indicates that we do not know 
the sign of the curvature term further implying that it is subdominant in 
Eq. © 

k 8ttG 

In an adiabatically expanding Universe, R ~ T^ 1 where T is the tempera- 
ture of the thermal photon background. Therefore the quantity 

k = -^-^r < T < 2 x 1CT 58 (12) 

R 2 T 2 3T 2 V ; 

is dimensionless and constant in the standard model. This is known as 
the curvature problem and can be resolved by a period of inflation. Before 
inflation, let us write R = Ri, T — Ti and R ~ T~ l . During inflation, 
R ~ T _1 ~ e Ht , where H is constant. After inflation, R = Rf 3> Ri but 
T = Tf = Tr < Ti where Tr is the temperature to which the Universe 
reheats. Thus R / T and k — > is not constant. But from Eqs. and 
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iJEJ if k -> then Q, -> 1, and since typical inflationary models contain 
much more expansion than is necessary, S7 becomes exponentially close to 
one. 

The inflationary prediction of SI = 1 is remarkably consistent with the 
CMB measurements discussed above. Furthermore, we know two things: 
Dark matter exists, since we don't see SI = 1 in luminous objects, and most 
(about 90%) of the dark matter is not baryonic. The latter conclusion is a 
result of our forthcoming discussion on BBN which constrains the baryon- 
to-photon ratio and hence SI b ■ Thus 1 — S7# is not only dark but also non- 
baryonic. Furthermore, the matter density is surely composed of several 
contributions: Sl m = Qb + Sl„ + S7 X where the latter represents the dark 
matter contribution. 

Another important piece of theoretical evidence for dark matter comes 
from the simple fact that we are living in a galaxy. The type of perturba- 
tions produced by inflation^! are, in most models, adiabatic perturbations 
{Sp/p tx ST/T), and I will restrict my attention to these. Indeed, the 
perturbations produced by inflation also hav e th e very nearly scale-free 
spectrum described by Harrison and Zeldovicrl "'. When produced, scale- 
free perturbations fall off as & cx I 2 (increase as the square of the wave 
number). At early times Sp/p grows as t until the time when the horizon 
scale (which is proportional to the age of the Universe) is comparable to 
I. At later times, the growth halts (the mass contained within the volume 
I 3 has become smaller than the Jean's mass) and ^ = 5 (roughly) inde- 
pendent of the scale I. When the Universe becomes matter dominated, the 
Jean's mass drops dramatically and growth continues as — cx R ~ 1/T. 

The transition to matter dominance is determined by setting the energy 
densities in radiation (photons and any massless neutrinos) equal to the 
energy density in matter (baryons and any dark matter) . For three massless 
neutrinos and baryons (no dark matter), matter dominance begins at 

T m = 0.22m B r? (13) 

and for r\ < 7 x 1CU 10 , this corresponds to T m < 0.14 eV. 

Because we are considering adiabatic perturbations, there will be 
anisotropies produced in the microwave background radiation on the order 
of ST/T r~j 5. The value of 8, the amplitude of the density fluctuations at 
horizon crossing, has been determined by COBE^J d = (5.7 ± 0.4) x 10~ 6 . 
Without the existence of dark matter, Sp/p in baryons could then achieve 
a maximum value of only Sp/p ~ A\S(T m /T ) < 2 x 10~ 3 A\, where 
T Q = 2.35 x 10 -4 eV is the present temperature of the microwave back- 
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ground and A\ ~ 1 — 10 is a scale dependent growth factor. The overall 
growth in 5p/p is too small to argue that growth has entered a nonlinear 
regime needed to explain the large value (10 5 ) of 5p/p in galaxies. 

Dark matter easily remedies this dilemma in the following way. The 
transition to matter dominance is determined by setting equal to each other 
the energy densities in radiation (photons and any massless neutrinos) and 
matter (baryons and any dark matter). While if we suppose that there 
exists dark matter with an abundance Y x — n x jn 1 (the ratio of the number 
density of x' s to photons) then 

T m = 0.22m x Y x (14) 

Since we can write m x Y x /GeV = fl x h 2 /(A x 10 7 ), we have T m /T Q — 2.4 x 
10 4 tt x h 2 which is to be compared with 600 in the case of baryons alone. 
The baryons, although still bound to the radiation until decoupling, now 
see deep potential wells formed by the dark matter perturbations to fall 
into and are no longer required to grow at the rate 5p/p oc R. 

With regard to dark matter and galaxy formation, all forms of dark 
matter are n ot e qual. They can be distinguished by their relative temper- 
ature at T m E2 Particles which are still largely relativistic at TL (like 
neutrinos or other particles with m x < 100 eV) have the property-^ that 
(due to free streaming) they erase perturbations out to very large scales 
given by the Jean's mass 

Mj = 3 x 10 18 ^- (15) 

Thus, very large scale structures form first and galaxies are expected to 
fragment out later. Particles with this property are termed hot dark matter 
particles. Cold particles (m x > 1 MeV) have the opposite behavior. Small 
scale structure forms first aggregating to form larger structures later. It is 
now well known that pure HDM cosmologies can not reproduce the observed 
large scale structure of the Universe. In contrast, CDM does much better. 
Current attention is focused on so-called ACDM cosmologies based on the 
fl\ — Cl m contribution to the curvature discussed above. 



2. Lecture 2 

2.1. Big Bang Nucleosynthesis 

The standard modeP^l f big bang nucleosynthesis (BBN) is based on the 
relatively simple idea of including an extended nuclear network into a ho- 
mogeneous and isotropic cosmology. Apart from the input nuclear cross 
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sections, the theory contains only a single parameter, namely the baryon- 
to-photon ratio, rj. Other factors, such as the uncertainties in reaction 
rates, and the neutron mean-li fe can be treated by standard statistical and 
Monte Carlo techniqueiP^EID The theory then allows one to make pre- 
dictions (with well-defined uncertainties) of the abundances of the light 
elements, D, 3 He, 4 He, and 7 Li. 



2.1.1. Theory 

Conditions for the synthesis of the light elements were attained in the early 
Universe at temperatures T > 1 MeV. In the early Universe, the energy 
density was dominated by radiation with 

P=^ + l + l^ (16) 

from the contributions of photons, electrons and positrons, and N v neutrino 
flavors (at higher temperatures, other particle degrees of freedom should be 
included as well). At these temperatures, weak interaction rates were in 
equilibrium. In particular, the processes 

n + e + <-> p + v e 
n + v e <-> p + e~ 

n «-> p + e~ + v e (17) 

fix the ratio of number densities of neutrons to protons. At T ^> 1 MeV, 
(n/p) ~ 1. 

The weak interactions do not remain in equilibrium at lower temper- 
atures. Freeze-out occurs when the weak interaction rate, ~ G^T 5 
falls below the expansion rate which is given by the Hubble parameter, 
H ~ VGWp - T 2 /M P , where M P = 1/v^jv - 1.2 x 10 19 GeV. The 
/^-interactions in eq. (|17|l freeze-out at about 0.8 MeV. As the tempera- 
ture falls and approaches the point where the weak interaction rates are 
no longer fast enough to maintain equilibrium, the neutron to proton ratio 
is given approximately by the Boltzmann factor, (n/p) ~ e ~ Am / T ^ 1/6, 
where Am is the neutron-proton mass difference. After freeze-out, free 
neutron decays drop the ratio slightly to about 1/7 before nucleosynthesis 
begins. 

The nucleosynthesis chain begins with the formation of deuterium by 
the process, p + n — > D + 7. However, because of the large number of 
photons relative to nucleons, rf 1 = n^/ne ~ 10 10 , deuterium production 
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is delayed past the point where the temperature has fallen below the deu- 
terium binding energy, Eb = 2.2 MeV (the average photon energy in a 
blackbody is E 1 ~ 2. .IT). This is because there are many photons in the 
exponential tail of the photon energy distribution with energies E > Eb 
despite the fact that the temperature or E 1 is less than Eb- The degree 
to which deuterium production is delayed can be found by comparing the 
qualitative expressions for the deuterium production and destruction rates, 

T p sa nscrv (18) 
T d « n 1 ave~ EB/T 

When the quantity f? _1 exp(— Eb/T) ~ 1, the rate for deuterium destruc- 
tion (D + 7 — * p + n) finally falls below the deuterium production rate and 
the nuclear chain begins at a temperature T ~ O.lMeV. 

The dominant product of big bang nucleosynthesis is He and its abun- 
dance is very sensitive to the (n/p) ratio 



Y P = « 0.25 (19) 

[1 + {n/p)} 



i.e., an abundance of close to 25% by mass. Lesser amounts of the other 
light elements are produced: D and 3 He at the level of about 1CP 5 by 
number, and 7 Li at the level of 10 -10 by number. 

The resulting abundances of the light element are shown in Figure||3 
over the range in 7710 = 10 10 ?7 between 1 and 10. The left plot shows the 
abundance of 4 He by mass, Y, and the abundances of the other three iso- 
topes by number. The curves indicate the central predictions from BBN, 
while the bands correspond to the uncertainty in the predicted abundances 
based primarily the un cert ainty in the input nuclear reactions as computed 
by Monte Carlo in ref. This theoretical uncertainty is shown explicitly 
in the right panel as a function of 7710 . The dark shaded boxes correspond 
to the observed abundances of 4 He and 7 Li and will be discussed below. 
The dashed boxes correspond to the ranges of the elements consistent with 
the systematic uncertainties in the observations. The broad band shows a 
liberal range for rjio consistent with the observations. At present, there is 
a general concordance between the theoretical predictions and the observa- 
tional data. 
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The light element abundances from big bang nucleosynthesis as a function of 



2.1.2. Abundances 

In addition to it BBN production, 4 He is made in stars, and thus co- 
produced with heavy elements. Hence the best sites for determining the 
primordial 4 He abundance are in metal-poor regions of hot, ionized gas in 
nearby external galaxies (extragalactic HII regions). Helium indeed shows a 
linear correlation with metallicity in these systems, and the extrapolatio n to 
zero metallicity gives the primordial abundance (baryonic mass fractionj^ 



Y n = 0.238 ±0.002 ±0.005. 



(20) 



Here, the first error is statistical and reflects the large sample of systems, 
whilst the second error is systematic and dominates. 

The systematic unc erta inties in these observations have not been thor- 
oughly explored to date^Sl. In particular, there may be reason to suspect 
that the above primordial abundance will be increased due to effects such 
as underlying stellar absorption in the HII regions. We note that other 
analyses give similar results: Y p = 0.244 ± 0.002 ± 0.005 49 and 0.239 
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±0.002 E@ 

The primordial 7 Li abundance comes from measurements in the atmo- 
spheres of primitive (Population II) stars in the stellar halo of our Galaxy. 
The 7 Li/H abundance is found to be constant for stars with low metallicity, 
indicating a primordial component, and a recent determination gives 

= (1.23 ± 0.06±g-g) x 10- 10 (95% CL), (21) 

where the sm all statistical error is overshadowed by systematic 
uncertainties^. The range (|21(l may, however, be underestimated, as a 
recent determination^! uses a different procedure to determine stellar at- 
mosphere parameters, and gives 7 Li/H p = (2.19 ± 0.28) x 10~ 10 . At this 
stage, it is not possible to determine which method of analysis is more 
accurate, indicating the likelihood that the upper systematic uncertainty 
in (|2*T]l has been underestimated. Thus, in order to obtain a conservative 
bound from 7 Li, we take the lower bound from l|21|) and the upper bound 
frorrP, giving 

7 J j 

9.0xl0" n <— <2.8xl0~ 10 . (22) 
H p 

Deuterium is measured in high-redshift QSO absorption line systems 
via its isotopic shift from hydrogen. In several absorbers of moderate col- 
umn density (Lyman-limit systems), D has been observed in multiple Ly- 
man transition ;] I I Restricting our attention to the three most reliable 
regionJ^, we find a weighted mean of 

^ = (2.9 ±0.3) x 10" 5 . (23) 
Hp 

It should be noted, however, that the x 2 per degree of freedom is rather poor 
( ~ 3.4), and that the unweighted dispersion of these data is ~ 0.6 x 10~ 5 . 
This already points to the dominance of systematic effects. Observation 
of D in systems with higher column density (damped systems) hnd lower 
D /hG22 at a level inconsistent with 123(1 , further suggesting that systematic 
effects dominate the error budget^. If all five available observations are 
used, we would find D/H = (2.6 ± 0.3) x 10 -5 with an even worse x 2 P er 
degree of freedom (~ 4.3) and an unweighted dispersion of 0.8. 

Because there are no known astrophysical sites for the production of 
deuterium, all observed D is assumed to be primordial^! As a result, any 
firm determination of a deuterium abundance establishes an u ppe r bound 
on rj which is robust. Thus, the recent measurements of D /H 1^ at least 
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provide a lower bound on D/H, D/H > 2.1 x 10 _5 (2er) and hence provide 
an upper bound to 77, 7710 < 7.3 and fish 2 < 0.027. 

Hclium-3 can be measured through its hypcrfinc emission in the radio 
band, and has been observed in HII regions in our Galaxy. These observa- 
tions fincP^Sl that there are no obvious trends in 3 He with metallicity and 
location in the Galaxy. There is, however, considerable scatter in the data 
by a factor ~ 2, some of which may be real. Unfortunately, the stellar and 
Galactic evolution of 3 He is not yet sufficiently well understood to confirm 
whether 3 He is increasing or decreasing from its primordial valud^l Con- 
sequently, it is unclear whether the observed 3 He abundance represents an 
upper or lower limit to the primordial value. Therefore, we can not use 3 He 
abundance as a constraint. 

By combining the predictions of BBN calculations with the abundances 
of D, 4 He, and 7 Li discussed above one can determine the the 95% CL 
region 4.9 < 7710 < 6.4, with the peak value occurring at 7710 = 5.6. This 
range corresponds to values of Ob between 

0.018 < n B h 2 < 0.023 (24) 

with a central value of fish 2 = 0.020. 

If we were to use only the deuterium abundance from Eq. 1231 one 
obtains the 95% CL range 5.3 < 7710 < 7.3, with the peak value occurring 
at 7710 = 5.9. This range corresponds to values of Qb between 

0.019 < n B h 2 < 0.027 (25) 

with a central value of Sl^/i 2 = 0.021. As one can see from a comparison 
with Table 1, these values are in excellent agreement with determinations 
from the CMB. 

2.2. Candidates 
2.2.1. Baryons 

Accepting the dark matter hypothesis, the first choice for a candidate 
should be something we know to exist, baryons. Though baryonic dark 
matter can not be the whole story if Cl m > 0.1, the identity of the dark 
matter in galactic halos, which appear to contribute at the level of fi ~ 0.05, 
remains an important question needing to be resolved. A baryon density of 
this magnitude is not excluded by nucleosynthesis. Indeed we know some 
of the baryons are dark since fl < 0.01 in the disk of the galaxy. 

It is interesting to note that until recently, there seemed to be some dif- 
ficulty in reconciling the baryon budget of the Universe. By counting the 
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visible contribution to f2 in stellar populations and the X-ray producing hot 
gas, Persic and Saluccf^ found only f2 v is — 0.003. A subsequent account- 
ing by Fukugita, Hogan and PeebleJSU found slightly more (fi ~ 0.02) by 
including the contribution from plasmas in groups and clusters. At high 
redshift on the other hand, all of the baryons can be accounted for. The 
observed opacity of the Ly a forest in QSO absorption spectra requires a 
large baryon density consistent with the determinations by the CMB and 

In galactic halos, however, it is quite difficult to hide large amounts of 
baryonic matter^" 2 * Sites for halo baryons that have been discussed include 
Hydrogen (frozen, cold or hot gas), low mass stars/ Jupiters, remnants of 
massive stars such as white dwarfs, neutron stars or black holes. In almost 
every case, a serious theoretical or observational problem is encountered. 

2.2.1.1 Hydrogen 

A halo predominately made of hydrogen (with a primordial admixture 
of He) is perhaps the simplest possibility. Hydrogen may however be 
present in a condensed snow-ball like state or in the form of gas. Aside 
from the obvious question of how do these snowballs get made, it is possible 
to show that their existence today requires them to be so large as to be 
gravitationally bounds 

Assuming that these objects are electrostatically bound, the average 
density of solid hydrogen is p s = 0.07 g cm -3 and the binding energy per 
molecule is about 1 eV. Given that the average velocity of a snowball is 
v ~ 250 kms" 1 , corresponding to a kinetic energy Ek ~ 600 eV, snowballs 
must be collisionless in order to survive. Requiring that the collision rate 
r c = n s av be less than t u ~ Y (t u is the age of the halo ~ age of the Universe) 
with n s = pn/rris and a = irr s 2 one finds r s > 2 cm and m s > 1 g, assuming 
a halo density pn = 1-7 x 10 -26 g cm" 3 (corresponding to 10 12 M Q in a 
radius of 100 kpc). However, collisionless snowballs also require that their 
formation occur when the overall density p = pn- In this case, snowballs 
could not have formed later than a redshift (1 + z) = 3.5 or when the 
microwave background temperature was 9.5K. At this temperature, there is 
no equilibrium between the gaseous and condensed state and the snowballs 
would sublimate. For a snowball to survive, r s > 10 16 cm is required 
making this no longer an electrostatically bound object. 

If snowballs sublimate, then we can consider the possibility of a halo 
composed of cold hydrogen gas. Because the collapse time-scale for the 
halo (< 10 9 yrs) is much less than the age of the galaxy, the gas must be 
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in hydrostatic equilibrium. Combining the equation of state 



P{r) = (2p(r)/m p )kT 



(26) 



where m p is the proton mass, with 



dP(r)/dr = — GM{r)p{r) 



(27) 



one can solve for the equilibrium temperature 



Gm„M(r) „ 
4kr 



(28) 



This is hot gas. As discussed earlier, hot gas is observed through X-ray 
emission. It is easy to show that an entire halo of hot gas would conflict 
severely with observations. Cooling of course may occur, but at the expense 
of star formation. 

2.2.1.2 Jupiter-like objects 

A very popular candidate for baryonic dark matter is a very low mass 
star or JLO. These are objects with a mass m < m = 0.08M Q , the mass 
nece ssar y to commence nuclear burning. Presumably there is a minimum 
mas m > rrimin = (0.004 — O.OO7)Af based on fragmentation, but the 
exact value is very uncertain. The contribution of these objects to the dark 
matter in the halo depends on how much mass can one put between m m i n 
and m and depends on the initial mass function (IMF) in the halo. 

An IMF is the number of stars formed per unit volume per unit mass 
and can be parameterized as 



In this parametrization, the Salpeter mass function corresponds to x = 1.35. 
Because stars are not observed with m < m , some assumptions about the 
IMF for low masses must be made. An example of the observed^ IMF in 
the solar neighborhood is shown in Fig. 

It is possible to use infrared observation^^ to place a lower limit on the 
slope, x, of the IMF in the halo of galaxies by comparing a mass-to-light 
ratio defined by 



where the total mass density in JLO's and low mass stars is given by 



(f> = Am-^ 1+X '> 



(29) 



Q = {pm/pL) Lq/Mq 



(30) 




(31) 
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Figure 7. The IMF in the solar neighborhood^! 



where ttiq = 0.75A/© is the mass of a giant. The luminosity density given 
by such a stellar distribution is 



PL = 



Lcpdm + pa 



(32) 



where L(m) is the luminosity of a star of mass m. pc is t he contribution 
to the luminosity densit y d ue to giant stars. The observed^ lower limits 
on Q translate to a limrs^ on x 



x > 1.7 



(33) 



with a weak dependence of m m i n . 

Unfortunately, one can not use Eq. I|33|) to exclude JLO's since we do 
not observe an IMF in the halo and it may be different from that in the 
disk. One can however make a comparison with existing observations, none 
of which show such a steep slope at low masses. Indeed, most observations 
leading to a determination of the IMF (such as the one shown in Fig. 0) 
show a turn over (or negative slope). To fully answer the questions regard- 
ing JLO's in the halo, one needs a better understanding of star formation 
and the IMF. For now, postulating the existence of a large fraction of JLO's 
in halo is rather ad-hoc. 
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Despite the theoretical arguments against them, JLO's or massive com- 
pact halo objects (MACHOs) are candidates which are testable by the grav- 
itational microlensing of stars in a neighboring galaxy such as the LMd^H. 
By observing millions of stars and examining their intensity as a function 
of time, it is possible to determine the presence of dark objects in our halo. 
It is expected that during a lensing event, a star in the LMC will have its 
intensity rise in an achromatic fashion over a period St ~ 3 \J M/.001M© 
days. Indeed, microlensing candidates have been found®. For low mass 
objects, those with M < 0.1M©, it appears however that the halo fraction 
of MACHOs is very small. The relative amount of machos in the halo is 
typically expressed in terms of an optical depth. A halo consisting 100% 
of machos would have an optical depth of r ~ 5 x 10~ 7 . The most recent 
results of the MACHO collaboratioiP^ indicate that r = 12t\ x 10~ 8 , cor- 
responding to a macho halo fraction of about 20% with a 95% CL range of 
8 - 50% based on 13 - 17 events. They also exclude a 100% macho halo at 
the 95% CL. The typical macho mass falls in the range 0.15 - 0.9 M©. The 
EROS collaboration has set even stronger limits having observed 5 events 
toward the LMC and 4 toward the SMO^f The observed optical depth 
from EROS1 is r = 4±f x 10~ 8 and from EROS 2 r = 6±g x 10~ 8 . They 
have excluded low mass objects (M < 0.1M©) to make up less than 10% 
of the halo and objects with 2 x 10~ 7 Af© < M < 1M Q to be less than 25% 
of the halo at the 95% CL. 

2.2.1.3 Remnants of Massive Stars 

Next one should possibility that the halo is made up of the dead stellar 
remnants of stars whose initial masses were M > 1M©. Briefly, the problem 
which arises in this context is that since at least 40% of the stars initial 
mass is ejected, and most of this mass is in the form of heavy elements, a 
large population of these objects would contaminate the disk and prevent 
the existence of extremely low metallicity objects (Z ~ 10~ 5 ) which have 
been observed. Thus either dust (from ej ecta) o r dead remnants would be 
expected to produce too large a metallicit y^^ ! Clearly star formation is a 
very inefficient mechanism for producing dark matter. Many generations of 
stars would would be required to cycle through their lifetimes to continually 
trap more matter in remnants. 

This question has been studied in more detaiM By allowing a variable 
star formation rate, and allowing a great deal of flexibility in the IMF, a 
search for a consistent set of parameters so that the halo could be described 
primarily in terms of dead remnants (in this case white dwarfs) was per- 
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formed. While a consistent set of chemical evolution parameters can be 
found, there is no sensible theory to support this choice. In such a model 
however, the dark matter is in the form of white dwar fs a nd the remaining 
gas is is heavily contaminated by Carbon and Nitroge Though it is not 
excluded, it is hard to understand this corner of parameter space as being 
realistic. 

2.2.1.4 Black Holes 

There are several possibilities for black holes as the dark matter in halos. 
If the black holes are primordial which have presumably formed before 
nucleosynthesis, they should not be counted as baryonic dark matter and 
therefore do not enter into the present discussion. If the black holes are 
formed as the final stage of star's history and its formation was preceded 
by mass loss or a supernovae, then the previous discussion on remnants 
of massive stars applies here as well. However, it is also possible that the 
black halos were formed directly from very massive stars (m > 100M Q ?) 
through gravitational instability with no mass Though there are 

limits due to overheating the disli^ and stellar systemd^l. I n this case I 
know of no argument preventing a sufficiently large population of massive 
black holes as baryonic dark matter (Of course, now the IMF must have 
m min > 1OOM .) 



2.2.2. Neutrinos 

Light neutrinos (m < 30 eV) are a long-time standard when it comes to non- 
baryonic dark matteffHl. Light neutrinos produce structure on large scales, 
and the natural (minimal) scale for structure clustering is given in Eq 
Hence neutrinos offer the natural possibility for large scale structure 
including filaments and voids. Light neutrinos are, however, ruled out as a 
dominan t fo rm of dark matter because they produce too much large scale 
structure^! Because the smallest non-linear structures have mass scale Mj 
and the typical galactic mass scale is ~ 10 12 Mq, galaxies must fragment 
out of the larger pancake-like objects. The problem with such a scenario is 
that galaxies form lat cffi^^ {z < 1) whereas quasars and galaxies are seen 
out to redshifts z > 6. 

In the standard model, the absence of a right-handed neutrino state 
precludes the existence of a neutrino mass (unless one includes non- 
renormalizable lepton number violating interactions such HHLL). By 
adding a right-handed state z/r, it is possible to generate a Dirac mass for 
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the neutrino, m v = as is the case for the charged lepton masses, 

where h u is the neutrino Yukawa coupling constant, and v is the Higgs 
expectation value. It is also possible to generate a Majorana mass for the 
neutrino when in addition to the Dirac mass term, m v VBPhi a term Mvrvr 
is included. If M 3> m„, the see-saw mechanism produces two mass eigen- 
states given by m Vl ~ m^/M which is very light, and m„ 2 ~ M which is 
heavy. The state v\ is a potential hot dark matter candidate as v% is in 
general not stable. 

The simplicity of the standard big bang model allows one to compute 
in a straightforward manner the relic density of any stable particle if that 
particle was once in thermal equilibrium with the thermal radiation bath. 
At early times, neutrinos were kept in thermal equilibrium by their weak 
interactions with electrons and positrons. As we saw in the case of the 
/3-interaction used in BBN, one can estimate the thermally averaged low- 
energy weak interaction scattering cross section 

(av) ~ ff 4 T 2 /m 4 w (34) 

for T <C m\v ■ Recalling that the number density scales as n oc T 3 , we 
can compare the weak interaction rate T ~ n{<rv), with the expansion rate 
given by eqs. (jHJ with 

'=(£^IX>)f^ s f>T)r< (35) 

Neutrinos will be in equilibrium when T w k > H or 

T 3 > y/8w 3 N/90 m^r/Mp (36) 

where M P = G N 1/2 = 1.22 x 10 19 GeV is the Planck mass. For N = 43/4 
(accounting for photons, electrons, positrons and three neutrino flavors) we 
see that equilibrium is maintained at temperatures greater than 0(1) MeV 
(for a more accurate calculation see^^J). 

The decoupling scale of 0(1) MeV has an important consequence on the 
final relic density of massive neutrinos. Neutrinos more massive than 1 MeV 
will begin to annihilate prior to decoupling, and while in equilibrium, their 
number density will become exponentially suppressed. Lighter neutrinos 
decouple as radiation on the other hand, and hence do not experience the 
suppression due to annihilation. Therefore, the calculations of the number 
density of light (m„ < 1 MeV) and heavy [m v > 1 MeV) neutrinos differ 
substantially. 
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The number of density of light neutrinos with m v < 1 MeV can be 
expressed at late times as 

p v = m u Y u n^ (37) 

where Y v — n v jn^ is the density of i/'s relative to the density of photons, 
which today is 411 photons per cm 3 . It is easy to show that in an adi- 
abatically expanding universe Y v = 3/11. This suppression is a result of 
the e + e~ annihilation which occurs after neutrino decoupling and heats the 
photon bath relative to the neutrinos. In order to obtain an age of the Uni- 
verse, t > 12 Gyr, one requires that the matter component is constrained 

by 

Qh 2 < 0.3. (38) 

From this o ne f inds the strong constraint (upper bound) on Majorana neu- 
trino massed^ 



TOtot 



m u < 28eV. (39) 



where the sum runs over neutrino mass eigenstates. The limit for Dirac 
neutrinos depends on the interactions of the right-handed states (see dis- 
cussion below). Given the discussion of the CMB results in the previous 
section, one could make a case that the limit on Slti 2 should be reduced by 
a factor of 2, which would translate in to a limit of 14 eV on the sum of the 
light neutrino masses. As one can see, even very small neutrino masses of 
order 1 eV, may contribute substantially to the overall relic density. The 
limit (|39|l and the corresponding initial rise in Q> v h 2 as a function of m v is 
displayed in the Figure |51 (the low mass end with m u < 1 MeV). 

Combining the rapidly improving data on key cosmological parameters 
with the better statistics from large redshift surveys has made it possible 
to go a step forward along this path. It is now possible to set stringent 
limits on the light neutrino mass density fl„h 2 , and he nce on the neutrino 
mass based on the power spectrum of the Ly a foresl 2 ^, mtot < 5.5 eV, 
and the limit is even stronger if the total matter density, tt m is less than 
0.5. Adding additional ob servation constraints from the CMB and galaxy 
clusters drops this limrJ^ to 2.4 eV. This limit has recently been improved 
by the 2dF Galaxy redshifJSSl survey by comparing the derived power spec- 
trum of fluctuations with structure formation models. Focussing on the 
the presently favoured ACDM model, the neutrino mass bound becomes 
mtot < 1-8 eV for f2 m < 0.5. When even more constraints such as HST 
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Figure 8. Summary of the relic density of Dirac neutrinos (solid) including 

possible neutrino asymmetry of r}u = 5 X 10" 11 (dotted). 



Key project data, supernovae type la data, and BBN are included^ the 
limit can be pushed to m to t < 0.3 eV. 

The calculation of the relic density for neutrinos more massive than ~ 1 
MeV, is substantially more involved. The relic density is now determined 
by the freeze-out of neutrino annihilations which occur at T < m„, after 
annihilations have begun to seriously reduce their number densitjl^l The 
annihilation rate is given by 

2 

T ann = {av) ann n v ~ ^L{ mv Tfl 2 e- m "l T (40) 
m% 

where we have assumed, for example, that the annihilation cross section 
is dominated by vv — > // via Z-boson exchange a and (av) ann ~ m1/m A z . 
When the annihilation rate becomes slower than the expansion rate of the 
Universe the annihilations freeze out and the relative abundance of neu- 
trinos becomes fixed. For particles which annihilate through approximate 
weak scale interactions, this occurs when T ~ m x /20. The number density 



a While this is approximately true for Dirac neutrinos, the annihilation cross section of 
Majorana neutrinos is p-wave suppressed and is proportional of the final state fermion 
masses rather than m„. 
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of neutrinos is tracked by a Boltzmann-like equation, 

^ = -3^n-<at,>(n 2 -n§) (41) 

where uq is the equilibrium number density of neutralinos. By defining the 
quantity / = rt/T 3 , we can rewrite this equation in terms of x, as 

The solution to this equation at late times (small x) yields a constant value 
of /, so that n oc T 3 . 

Roughly, the solution to the Boltzmann equation goes as Y v ~ / ~ 
(m{av} a nn)~ 1 and hence Q, v h 2 ~ (crv) ann , so that parametrically <>,,//• 
1/m 2 ,. As a result, the constraint l|38l) now leads to a lower boundSQEIESl 
on the neutrino mass, of about m v > 3 — 7 GeV, depending on whether 
it is a Dirac or Majorana neutrino. This bound and the corresponding 
downward trend Q v h 2 ~ can again be seen in Figure El The result 

of a more detailed calculation is shown in Figure [5] for the case of a 
Dirac neutrino. The two curves show the slight sensitivity on the temper- 
ature scale associated with the quark-hadron transition. The result for a 
Majorana mass neutrino is qualitatively similar. Indeed, any particle with 
roughly weak scale cross-sections will tend to give an interesting value of 

nh 2 ~ 1. 

The deep drop in tt u h 2 , visible in Figure [5] at around m v = Mz/2, is 
due to a very strong annihilation cross section at Z-boson pole. For yet 
higher neutrino masses the Z-annihilation channel cross section drops as 
~ 1/m 2 , leading to a brief period of an increasing trend in f2„/i 2 . However, 
for TJijj ttiw the cross section regains its parametric form (o~v^ ann ^ m^, 
due to the opening up of a new annihilation channel to W^-boson paira^l 
and the density drops again as Vl u h 2 ~ 1/m 2 . The tree level M^-channel 
cross section breaks the unitarity at around O(few) Te Vjj^p however, and 
the full cross section must be bound by the unitarity limilPSI. This behaves 
again as 1/m 2 , whereby f2„ft. 2 h as to st art increasing again, until it becomes 
too large again at 200-400 TeV l ^5 | 94 | ^ or p erna p S somewhat earlier as the 
weak interactions become strong at the unitarity breaking scale). 

If neutrinos are Dirac particles, and have a nonzero asymmetry the relic 
density could be governed by the asymmetry rather than by the annihilation 
cross section. Indeed, it is easy to see that the neutrino mass density 
corresponding to the asymmetry rj u = (n v — n ij )/n 1 is given by^l 

p = m v r\ v n 1 , (43) 
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Figure 9. The relic density of heavy Dirac neutrinos due to annihilations*^. The curves 
are labeled by the assumed quark-hadron phase transition temperature in MeV. 



which implies 



where ?7„io = 10 10 ?7„. The behaviour of the energy density of neutrinos 
with an asymmetry is shown by the dotted line in the Figure |S] At low 
m„, the mass density is dominated by the symmetric, relic abundance of 
both neutrinos and antineutrinos which have already frozen out. At higher 
values of m„, the annihilations suppress the symmetric part of the relic 
density until Vl v h 2 eventually becomes dominated by the linearly increasing 
asymmetric contribution. In the figure, we have assumed an asymmetry of 
r\ v ~ 5 x 10 -11 for neutrinos with standard weak interaction strength. In 
this case, f2„/i 2 begins to rise when m„ > 20 GeV. Obviously, the bound 
(|3*5|l is saturated for m v = 75 GeV/^io- 

Based on the leptonic and invisible width of the Z boson, experi- 
ments at LEP have determined that the number of neutrinos is N v = 
2.9841 ± 0.0083 1*^1 Conversely, any new physics must fit within these 
brackets, and thus LEP excludes additional neutrinos (with standard weak 
interactions) with masses m v < 45 GeV. Combined with the limits dis- 
played in Figures |H1 and we see that the mass density of ordinary heavy 
neutrinos is bound to be very small, £l„h 2 < 0.001 for masses m v > 45 
GeV up to m v ~ 0(100) TeV. Lab constraints for Dirac neutrinos are 



n v h 2 



0.004^io (m„/GeV). 



(44) 
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available^, excluding neutrinos with masses between 10 GeV and 4.7 TeV. 
This is significant, since it precludes the possibility of neutrino dark matter 
based on an asymmetry between v and v * u . Majorana neutrinos are ex- 
cluded as dark matter since fl u h a 2 < 0.001 for m v > 45 GeV and are thus 
cosmologically uninteresting. 

A bound on neutrino masses even stronger than Eqn. (|39|) can be 
obtained from the recent observations of active-active mixing in both 
solar- and atmospheric neutrino experiments. The inferred evidence for 
Vjx — v T and v e — v^, T mixings are on the scales m 2 ~ 1 — 10 x 10 -5 and 
to 2 ~ 2 — 5 x 10 -3 . When combined with the upper bound on the electron- 
like neutrino mass m v < 2.8 eV 2S| and the LEP-limit on the number of 
neutrino species, one finds the constraint on the sum of neutrino masses: 

0.05 eV < mtot £ 8.4 eV. (45) 

Conversely, the experimental and observational data then implies that the 
cosmological energy density of all light, weakly interacting neutrinos can 
be restricted to the range 

0.0005 < n u h 2 < 0.09. (46) 

Interestingly there is now also a lower bound due to the fact that at least 
one of the neutrino masses has to be larger than the scale m 2 ~ 10~ 3 eV 2 
set by the atmospheric neutrino data. Combined with the results on relic 
mass density of neutrinos and the LEP limits, the bound (|46|) implies that 
the ordin ary weakly interacting neutrinos, once the standard dark matter 
candidate"*!, can be ruled out completely as a dominant component of the 
dark matter. 

If instead, we consider right-handed neutrinos, we have new possibilities. 
Right-handed interactions are necessarily weaker than standard left-handed 
interactions implying that right-handed neu trino s decouple early and today 
are at a reduced temperature relative to vl 

(Zk)3 = 43 (47) 
{ T 7 > AN(Td) K ' 

As such, for T dR > 1 MeV, n UR /n„ L = (T VR /T UL ) S < 1. Thus the abun- 
dance of right-handed neutrinos can be written as 

^ = ^ = (^)&) 3 «^ (48) 
n 7 11 1 VL 11 

In this case, the previous bound 139fl on neutrino masses is weakened. For a 
suitably large scale for the right-handed interactions, right-handed neutrino 
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masses may be as large as a few keV* '' '. Such neutrinos make excellent 
warm dark matter candi dates , albeit the viable mass range for galaxy for- 
mation is quite restrictecP^I. 

2.2.3. Axions 

Due to space limitations, the discussion of this candidate will be very brief. 
Axions are pseu do-Goldstone bosons which arise in solving the strong CP 
problenjnEIESI & global U(l) Peccei-Quinn symmetry. The invisible 
axio is associated with the flat direction of the spontaneously broken 
PQ symmetry. Because the PQ symmetry is also explicitly broken (the 
CP violating OFF coupling is not PQ invariant) the axion picks up a small 
mass similar to pion picking up a mass when chiral symmetry is broken. 
We can expect that m a ~ Tn^f w / f a where f a , the axion decay constant, is 
the vacuum expectation value of the PQ current and can be taken to be 
quite large. If we write the axion field as a = f a 9, near the minimum, the 
potential produced by QCD instanton effects looks like V ~ m 2 a 6 2 f%. The 
axion equations of motion lead to a relatively stable oscillating solution. 
The energy density stored in the oscillations exceeds the critical densitj^^ 
unless f a < 10 12 GeV. 

Axions may also be emitted stars and supernov d^^ l In supernovae, 
axions are produced via nucleon-nucleon bremsstrahlung with a coupling 
gAN oc mN I fa- As was noted above the cosmological density limit requires 
fa ^ 10 12 GeV. Axion emission from red giants imply-221 j a > \qW Q e y 
(though this limit depe nds o n an adjustable axion-electron coupling), the 
supernova limit requires^^ /„ > 2 x 10 11 GeV for naive quark model 
couplings of the axion to nucleons. Thus only a narrow window exists for 
the axion as a viable dark matter candidate. 

3. Lecture 3: Supersymmetric Dark Matter 

Although there are many reasons for considering supersymmetry as a can- 
didate exte nsion to the standard model of strong, weak and electromagnetic 
interactions^^, o ne of the most compelling is its role in understanding the 
hierarchy problerr l^^^ l namely, why/how is mw *C Mp. One might think 
naively that it would be sufficient to set mw ^ Mp by hand. However, 
radiative corrections tend to destroy this hierarchy For example, one-loop 
diagrams generate 




(49) 
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where A is a cut-off representing the appearance of new physics, and the 
inequality in (|49|l applies if A ~ 10 3 TeV, and even more so if A ~ incur ~ 
10 16 GeV or ~ Mp ~ 10 19 GeV. If the radiative corrections to a physical 
quantity are much larger than its measured values, obtaining the latter 
requires strong cancellations, which in general require fine tuning of the 
bare input parameters. However, the necessary cancellations are natural in 
supersymmetry, where one has equal numbers of bosons and fermions with 
equal couplings, so that Ij49(l is replaced by 



Sm 2 w =C?(^) \m%-m%\ . (50) 



The residual radiative correction is naturally small if \m 2 B — m F \ < 1 TeV 2 . 

In order to justify the absence of interactions which can be responsible 
for extremely rapid proton decay, it is common in the minimal supersym- 
metric standard model (MSSM) to assume the conservation of R-parity. If 
R-parity, which distinguishes between "normal" matter and the supersym- 
metric partners and can be defined in terms of baryon, lepton and spin as 
R — (— i) 3S + L + 2S ; i s unbroken, there is at least one supersymmetric par- 
ticle (the lightest supersymmetric particle or LSP) which must be stable. 
Thus, the minimal model contains the fewest number of new particles and 
interactions necessary to make a consistent theory. 

There are very strong constraints, however, forbidding the existen ce of 
stable or long lived particles which are not color and electrically neutrainH. 
Strong and electromagnetically interacting LSPs would become bound with 
normal matter forming anomalously heavy isotopes. Indeed, there are very 
strong uppe r limits on the abundances, relative to hydrogen, of nuclear 
isotopeJS^, n/n H < 10" 15 to 10" 29 for 1 GeV < m < 1 TeV. A strongly 
interacting stable relic is expected to have an abundance n/nn Ss 10~ 10 
with a higher abundance for charged particles. 

There are relatively few supersymmetric candidates which are not col- 
ored and are electrically neutral. The sneutrino x ^* 1 is one possibility, but 
in the MSSM, it has been excluded as a dark matter candidate by direc 
and indire clJlillsearches. In fact, one can set an accelerator based limit on 
the sneutrino mass from neutrino counting, mp > 44.7 GeV Hli3 i n this 
case, the direct relic searches in underground low-background experiments 
require mp > 20 TeV l^&l. Another possibility is the gravitino which is 
probably the most difficult to exclude. I will concentrate on the remaining 
possibility in the MSSM, namely the neutralinos. 
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3.1. Parameters 

The most general version of the MSSM, despite its minimality in particles 
and interactions contains well over a hundred new parameters. The study 
of such a model would be untenable were it not for some (well motivated) 
assumptions. These have to do with the parameters associated with super- 
symmetry breaking. It is often assumed that, at some unification scale, all 
of the gaugino masses receive a common mass, mx/2- The gaugino masses 
at the weak scale are determined by running a set of renormalization group 
equations. Similarly, one often assumes that all scalars receive a common 
mass, mo, at the GUT scale. These too are run down to the weak scale. 
The remaining parameters of importance involve the Higgs sector. There is 
the Higgs mixing mass parameter, and since there are two Higgs doublets 
in the MSSM, there are two vacuum expectation values. One combination 
of these is related to the Z mass, and therefore is not a free parameter, 
while the other combination, the ratio of the two vevs, tan (3, is free. 

If the supersymmetry breaking Higgs soft masses are also unified at the 
GUT scale (and take the common value mo), then [i and the physical Higgs 
masses at the weak scale are determined by electroweak vacuum conditions 
(/i is determined up to a sign). This scenario is often referred to as the 
constrained MSSM or CMSSM. Once these parameters are set, the entire 
spectrum of sparticle masses at the weak scale can be calculated 13 . In Fig. 
1101 an example of the running of the mass parameters in the CMSSM is 
shown. Here, we have chosen TOi/ 2 = 250 GeV, m = 100 GeV, tan/3 = 
3, Aq = 0, and (i < 0. Indeed, it is rather amazing that from so few 
input parameters, all of the masses of the supersymmetric particles can 
be determined. The characteristic features that one sees in the figure, are 
for example, that the colored sparticles are typically the heaviest in the 
spectrum. This is due to the large positive correction to the masses due 
to Qi3 in the RGE's. Also, one finds that the B (the partner of the U(1)y 
gauge boson), is typically the lightest sparticle. But most importantly, 
notice that one of the Higgs mass 2 , goes negative triggering electroweak 
symmetry breakin g^ ^ 1 (The negative sign in the figure refers to the sign 
of the mass 2 , even though it is the mass of the sparticles which are depicted.) 



b There are in fact, additional parameters: the supersymmetry-breaking tri-linear masses 
A (also assumed to be unified at the GUT scale) as well as two CP violating phases 9 M 
and 9a- 
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Figure 10. RG evolution of the mass parameters in the CMSSM. 

3.2. Neutralinos 

There are four neutralinos, each of which is a linear combination of the R = 
— 1 neutral fermionJUil the wino W 3 , the partner of the 3rd component 
of the SU(2)l gauge boson; the bino, B; and the two neutral Higgsinos, 
Hi and H.%. Assuming gaugino mass universality at the GUT scale, the 
identity and mass of the LSP are determined by the gaugino mass m\/2, (J-, 
and tan [3. In general, neutralinos can be expressed as a linear combination 

X = uB + (3W 3 +jHx + SH2 (51) 

The solution for the coefficients a, /?, 7 and S for neutralinos that make up 
the LSP can be found by diagonalizing the mass matrix 
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where Mi(M 2 ) is a soft supersymmetry breaking term giving mass to the 
U(l) (SU(2)) gaugino(s). In a unified theory Mi = M 2 — min at the 
unification scale (at the weak scale, Mi ~ |^-M 2 ). As one can see, the 
coefficients a, f3, 7, and 5 depend only on mi/2, fJ-, and tan/3. 
In Figure 

ED 117 . 

regions in the M%, fj, plane with tan/3 = 2 are shown 
in which the LSP is one of several nearly pure states, the photino, 7, the 
bino, B, a symmetric combination of the Higgsinos, -ff(i2), or the Higgsino, 
S — sin (3Hi +cos /3i?2- The dashed lines show the LSP mass contours. The 
cross hatched regions correspond to parameters giving a chargino (W , H^) 
state with mass < 45GeV and a s suc h are excluded by LEP " I This 
constraint has been extended by LEP^EH and is shown by the light shaded 
region and corresponds to regions where the chargino mass is < 103.5 GeV. 
The newer limit does not extend deep into the Higgsino region because 
of the degeneracy between the chargino and neutralino. Notice that the 
parameter space is dominated by the B or H12 pure states and that the 
photino only occupies a small fraction of the parameter space, as does 
the Higgsino combination S. Both of these light states are experimentally 
excluded. 

3.3. The Relic Density 

The relic abundance of LSP's is determined by solving the Boltzmann 
equation for the LSP number density in an expanding Universe. The 
techniquJ^ use d is similar to that for computing the relic abundance of 
massive neutrinos^! The relic density depends on additional parameters in 
the MSSM beyond mi/2,/i, and tan/3. These include the sfermion masses, 
and the Higgs pseudo-scalar mass, tua, derived from mo (and mj^)- To 
determine the relic density it is necessary to obtain the general annihilation 
cross-section for neutralinos. In much of the parameter space of interest, 
the LSP is a bino and the annihilation proceeds mainly through sfermion 
exchange. Because of the p-wave suppression associated with Majorana 
fermions, the s-wave part of the annihilation cross-section is suppressed by 
the outgoing fermion masses. This means that it is necessary to expand the 
cross-section to include p-wave corrections which can be expressed as a term 
proportional to the temperature if neutralinos are in equilibrium. Unless 
the neutralino mass happens to lie near near a pole, such as m x ~ mz/2 
or TOfr/2, in which case there are large contributions to the annihilation 
through direct s-channel resonance exchange, the dominant contribution 
to the BB annihilation cross section comes from crossed t-channel sfermion 
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Figure 11. Mass contours and composition of nearly pure LSP states in the MSSMliiH 



exchange. 

Annihilations in the early Universe continue until the annihilation rate 
r ~ avn x drops below the expansion rate. The calculation of the neutralino 
relic density proceeds in much the same way as discussed above for neutrinos 
with the appropriate substitution of the cross section. The final neutralino 
relic dens ity e xpressed as a fraction of the critical energy density can be 
written a a ' 



VL x h l 



1.9 x 10 



-li 



tJ f 



GeV 



ax f 



\bx) 



(53) 



where (T x /T 7 ) 3 accounts for the subsequent reheating of the photon tem- 
perature wi th re spect to x, due to the annihilations of particles with mass 
m < x/m x ^^. The subscript / refers to values at freeze-out, i.e., when 
annihilations cease. The coefficients a and b are related to the partial wave 
expansion of the cross-section, av = a + bx + . . .. Eq. I)53h results in a very 
good approximation to the relic density expect near s-channel annihilation 
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poles, thresholds and in regions where the LSP is nearly degenerate with 
the next lightest supersymmetric participial 

When there are several particle species i, which are nearly degenerate 
in mass, co-annihilations are important. In this cas the rate equation 
1|41JI still applies, provided n is interpreted as the total number density, 

n = J2 n i > ( 54 ) 



n n as the total equilibrium number density. 

E 



n = y~] n ,i , (55) 



and the effective annihilation cross section as 



(CTcffWrcl/ = 



—2 (o-y«rel) ■ (56) 



y u 



In eq. Q42[l. m x is now understood to be the mass of the lightest sparticle 
under consideration. 

Note that this implies that the ratio of relic densities computed with 
and without coannihilations is, roughly, 

( . 7) 




fi V &o J \ x 



where a = a + bx/2 and sub- and superscripts denote quantities com- 
puted ignoring coannihilations. The ratio Xf/xf ~ 1 + Xj la(g e sa e s / g\ao) , 
where g c g = ^ gi (wj/m i) 3 / 2 e~^ mi ~ mi ^ T . For the case of three degener- 
ate slepton NLSPsQZD g cS = J2i9* = 8 and x °f/ x f ~ I- 2 - The effects of 
co-annihilations are discussed below. 

3.4. Phenomenological and Cosmological Constraints 

For the cosmological limits on the relic density I will assume 

0.1 < Q x h 2 < 0.3. (58) 

The upper limit being a conservative bound based only on the lower limit 
to the age of the Universe of 12 Gyr. Indeed, most analyses indicate that 
^matter 0.4 — 0.5 and thus it is very likely that £l x h 2 < 0.2 (cf. the 
CMB results in Table 1). One should note that smaller values of Vt x h 2 are 
allowed, since it is quite possible that some of the cold dark matter might 
not consist of LSPs. 
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The calculated relic density is found to have a relevant cosmological 
density over a wide range of susy parameters. For all values of tan /?, 
there is a 'bulk' region with relatively low values of mi/a and mo where 
0.1 < Vt x h? < 0.3. However there are a number of regions at large values 
of mi/2 and/or mo where the relic density is still compatible with the 
cosmological constraints. At large values of mi/a, the lighter stau, becomes 
nearly degenerate with the neutralino and co-annihilations between these 
particles must be taken into account! 1 ' For non-zero value s of An , 
there are new regions for which \ — t coannihilations are importan t At 
large tan/?, as one increases m^, the pseudo-scalar mass, rriA begins to 
drop so that there is a wide funnel-like region (at all value s of mn) s uch that 
2m x rriA and s-channel annihilations become important^^ESl Finally, 
there is a region at very high mo where the value of /i begins to fall and 
the L SP b ecomes more Higgsino-like. This is known as the 'focus point' 
region 120 . 

As an aid to the assessment of the prospects for detecting sparticles at 
different accelerators, benchmark sets of supersymmetric parameters have 
often been found useful, since they provide a focus for concentrated discus- 
sion. A set of proposed post-LEP benchmark scenario in the CMSSM 
are illustrated schematically in Fig. ^] Five of the chosen points are in 
the 'bulk' region at small m^ 2 and mo, four are spread along the coan- 
nihilation 'tail' at larger m^ 2 for various values of tan/3. This tail runs 
along the shaded region in the lower right corner where the stau is the LSP 
and is therefore excluded by the constraints against charged dark matter. 
Two points are in rapid-annihilation 'funnels' at large m^ 2 and mo. At 
large values of mo, the focus-point region runs along the boundary where 
electroweak symmetry no longer occurs (shown in Fig. 1121 as the shaded 
region in the upper left corner). Two points were chosen in the focus-point 
region at large mo. The proposed points range over the allowed values of 
tan (3 between 5 and 50. The light shaded region corresponds to the portion 
of parameter space where the relic density f2 x /i 2 is between 0.1 and 0.3. 

The effect of coannihilations is to create an allowed band about 25-50 
GeV wide in mo for ^ 1400 GeV, which tracks above the m^ = m x 

contour. Along the line = m x , R w 10, from (|57|) 121 . As mo increases, 
the mass difference increases and the slepton contribution to <r e ff falls, and 
the relic density rises abruptly. This effect is seen in Fig. ^| The light 
shaded region corresponds to 0.1 < tth 2 < 0.3. The dark shaded region 
has m^ < m x and is excluded. The light dashed contours indicate the 
corresponding region in f2/i 2 if one ignores the effect of coannihilations. 



37 



e 

s 




Figure 12. Schematic overview of the CMSSM benchmark points proposed in^^3 The 
points are intended to illustrate the range of available possibilities. The labels correspond 
to the approximate positions of the benchmark points in the (mi^^mo) plane. They 
also span values of tan /3 from 5 to 50 and include points with fj, < 0. 

Neglecting coannihilations, one would find an upper bound of ~ 450 GeV 
on toi/2, corresponding to an upper bound of roughly 200 GeV on rrig. 
Instead, values of m 1 / 2 up to ~ 1400 GeV are allowed corresponding to an 
upper bound of ~ 600 GeV onmj. 

The most important phenomenological constraints are also shown 
schematically in Figure 1121 These include the constra int p rovided by the 
LEP lower limit on the Higgs mass: m# > 114.4 GeV I * \ This holds in 
the Standard Model, for the lightest Higgs boson h in the general MSSM 
for tan/3 < 8, and almost always in the CMSSM for all tan/3. Since mh 
is sensitive to sparticle masses, particularly my, via loop corrections, the 
Higgs limit also imposes important constraints on the CMSSM parameters, 
principally as seen by the dashed curve in Fig. 1121 



38 




Figure 13. The light-shaded 'bulk' area is the cosmologically preferred region with 
0.1 < £lh? < 0.3. The light dashed lines show the location of the cosmologically pre- 
ferred region if one ignores coannihilations with the light sleptons. In the dark shaded 
region in the bottom right, the LSP is the fi, leading to an unacceptable abundance 
of charged dark matter. Also shown is the isomass contour rn ± = 104 GeV and 
m/, = 110, 114 GeV, as well as an indication of the slepton bound from LEP. 



The constraint imposed by measurements of b — > sj ' also exclude 
small values of vriyii. These measurements agree with the Standard Model, 
and therefore provide bounds on MSSM particles, such as the chargino and 
charged Higgs masses, in particular. Typically, the b — > constraint is 
more important for /i < 0, but it is also relevant for fi > 0, particularly when 
tan/3 is large. The BNL E821 experiment reported last year a new mea- 
surement of a M = 5 (ffyi — 2) which deviated by 2.6 standard deviations from 
the best Standard Model prediction available at that tim J^^^ i However, it 
had been realized that the sign of the most importan t pse udoscalar-meson 
pole part of the light-by-light scattering contributio n ^ • to the Standard 
Model prediction should be reversed, which reduces the apparent experi- 
mental discrepancy to about 1.6 standard deviations (<5a M x 10 10 = 26± 16). 
The largest contribution to the errors in the comparison with theory was 
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thought to be the statistical erro r of the experiment, which has been signif- 
icantly reduced just recently^IS The world average of a M = — 2) now 
deviates by (33.9 ± 11.2) x 1CP 10 from the Standard Model calculation of 
Davier et al! fusing e + e~ data, and by (1 7 ± 11) x 10~ 10 from the Stan- 
dard Model calculation of Davier et al! based on r decay data. Other 
recent analyses of the e + e~ data yield similar results. On the subsequent 
plots, the formal 2-a range 11.5 x 10~ 10 < Sa^ < 56.3 x 10~ 10 is displayed. 
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Figure 14. Compilation of phenomenological constraints on t he CMSSM for tan /3 = 
10, At > 0, assuming Ao = 0, m t = 175 GeV and m b (m b )^j = 4.25 GeV . The near- 
vertical lines are the LEP limits "i x ± = 103.5 GeV (dashed and black^^l and = 
114.1 GeV ( dotted and red ^^ l Also, in the lower left corner we show the mg = 99 GeV 
contour In the dark (brick red) shaded regions, the LSP is the charged ri, so this 
region is excluded. The light (turquoise) shaded areas are the cosmologically preferred 
regions wi th 0.1 < Qh 2 < 0.3 l I The medium (dark green) shaded regions are excluded 
by b -> 57 11221 The shaded (pink) region in the upper right delineates the 2 a range of 
g M — 2. The dashed curves within this region correspond to the 1 — cr bounds. 



10 



Following a previous analysii^^IES] j n Fignre H^l the TO1/2 — mo param- 
eter space is shown for tan [3 — 10. The dark shaded region (in the lower 
right) corresponds to the parameters where the LSP is not a neutralino but 
rather a tr. The cosmologically interesting region at the left of the figure 
is due to the appearance of pole effects. There, the LSP can annihilate 
through s-channel Z and h (the light Higgs) exchange, thereby allowing a 
very large value of mo. However, this region is excluded by phenomeno- 
logical constraints. Here one can see clearly the coannihilation tail which 
extends towards large values of m 1 / 2 - In addition to the phenomenological 
constraints discussed above, Figure IT4l also shows the current experimental 
constraints on the CMSSM parameter space due t o the limit m x ± > 103.5 
GeV provided by chargino searches at LEP LEP has also provided 

lowe r limits on slepton masses, of which the strongest is mg > 99 GeV 
E23 This is shown by dot-dashed curve in the lo wer l eft corner of Fig. ^] 
Similar results have been found by other analyseJi^l 

As one can see, one of the most important phenomenological constraint 
at this value of tan /3 is due to the Higgs mass (shown by the nearly verti- 
cal dot-da shed curve). The theoretical Higgs masses were evaluated using 
FeynHiggJEHI which is estimated to have a residual uncertainty of a couple 
of GeV in m/j. The region excluded by the b — > sj constraint is the dark 
sh aded (green) region to the left of the plot. As many authors have pointed 
on a discrepancy between theory and the BNL experiment could well 
be explained by supersymmetry. As seen in Fig. 1141 this is particularly easy 
if fi > 0. The medium (pink) shaded region in the figure corresponds to 
the overall allowed region (2a) by the new experimental result. 

As discussed above, another mechanism for extending the allowed 
CMSSM region to large m x is rapid annihilation via a direct-channel pole 
when m x ~ \"tnA | 1^4 | 125 | This m ay yield a 'funnel' extending to large 
7774/2 and mo at large tan/3, as seen in Fig. 1151 

In principle, the true input parameters in the CMSSM are: /i, 777 1; 777 2 , 
and B, where mi and 777,2 are the Higgs soft masses (in the CMSSM 
777 1 = m 2 = 77io and B is the susy breaking bilinear mass term). In this 
case, the electroweak symmetry breaking conditions lead to a prediction of 
Mz,tan/3 ,and itia- Since we are not really interested in predicting Mz, 
it is more useful to assume instead the following CMSSM input param- 
eters: M£,777i,m 2 , and tan/3 again with mi = 7?7 2 = mo. In this case, 
one predicts /J.,B, and itia- However, one can generalize the CMSSM case 
to include non-universal Higgs masse^^EO (NUHM), in which case the 
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Figure 15. As in Fig. El for tan = 50. 



input parameters become:M^, /x, ttia, and tan/3 and one predicts m\,rri2, 
and B. 

The NUHM parameter space was recently analyze and a sample of 
the results found is shown in Fig. El While much of the cosmologically 
preferred area with fi < is excluded, there is a significant enhancement in 
the allowed parameter space for /i > 0. 



3.5. Detection 

Because the LSP as dark matter is present locally, there are many avenues 
for pursuing dark matter detection. Direct detection techniques rely on 
an ample neutralino-nucleon scattering cross-section. The effective four- 
fermion lagrangian can be written as 
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Figure 16. Compilations of phenomenological constraints on the MSSM with NUHM 
in the (p, itia) plane for tan /3 = 10 and mo = 100 GeV, Bljn = 300 GeV, assuming 

Aq = 0, mt = 175 GeV and mj,(mj,)c^| = 4.25 GeV. The shading is as described in 
Fig. 1141 The (blue) solid line is the contour m x = near which rapid direct-channel 

annihilation suppresses the relic density. The dark (black) dot-dashed line indicates when 
one or another Higgs mass-squared becomes negative at the GUT scale: only lower |^t| 
and larger tua values are allowed. The crosses denote the values of fi and rriA found in 
the CMSSM. 

+ a 5 txxqti 5 qt + a6ixi b xmi (59) 

However, the terms involving an, an, a^i, and a§i lead to velocity depen- 
dent elastic cross sections. The remaining terms are: the spin dependent 
coefficient, a%i and the scalar coefficient a^i. Contributions to a2i are 
predominantly through light squark exchange. This is the dominant chan- 
nel for binos. Scattering also occurs through Z exchange but this channel 
requires a strong Higgsino component. Contributions to a%i are also dom- 
inated by light squark exchange but Higgs exchange is non-negligible in 
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most cases. 
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Figure 17. Spin-independent cross sections in the (m-i/2,rno) planes for (a) tan/3 = 
10, fj, < 0, (b) tan/3 = 10, fi > 0. The double dot-dashed (orange) curves are contours of 
the spin-independent cross section, differing by factors of 10 (bolder) and interpolating 
factors of 3 (finer - when shown). For example, in (b), the curves to the right of the 
one marked 10 — 9 pb correspond to 3 X 10~ 10 pb and 10~ 10 pb. 



Fig. 1171 displays contours of the spin-independent cross section for the 
elastic scattering of the LSP x on protons in the m^jWo planes for (a) 
tan/3 = 10, /it < 0, (b) tan/3 = 10, a* > E2. The double dot-dashed 
(orange) lines are contours of the spin-independent cross section, and the 
contours asi = 10~ 9 pb in panel (a) and agi — 10~ 12 pb in panel (b) are 
indicated. The LEP lower limits on rrih and m x ± , as well as the experi- 
mental measurement of b — > sj for ft < 0, tend to bound the cross sections 
from above, as discussed in more detail below. Generally speaking, the 
spin-independent cross section is relatively large in the 'bulk' region, but 
falls off in the coannihilation 'tail'. Also, we note also that ther e is a stro ng 
cancellation in the spin-independent cross section when /i < 1 142 | 143 ag 
seen along strips in panel (a) of Fig. 1171 where m 1 / 2 ~ 500 GeV. In the 
cancellation region, the cross section drops lower than 10 -14 pb. All these 
possibilities for suppressed spin-independent cross sections are disfavoured 
by the data on — 2, which favour values of mi/2 and mo that are not 
very large, as well as \x > 0, as seen in panel (b) of Fig. El Thus — 2 
tends to provide a lower bound on the spin-independent cross section. 

Fig. lTM a) illustrates the effect on the cross sections of each of the princi- 
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Figure 18. Allowed ranges of the cross sections for tan/3 = 10 (a) fi > for spin- 
independent elastic scattering. The solid (blue) lines indicate the relic density constraint, 
the dashed (black) lines the b — > sj constraint, the dot-dashed (green) lines the mj, 
constraint, and the dotted (red) lines the g^ — 2 constraint. The shaded (pale blue) 
region is allowed by all the constraints, (b) The allowed ranges of the spin-independent 
cross section for fi > 0. The darker solid (black) lines show the upper limits on the 
cross sections obtained from and b — » S7, and (where applicable) the lighter solid 
(red) lines show the lower limits suggested by g M — 2 and the dotted (green) lines the 
lower limits from the relic density. 



pal phenomenological constraints, for the particular case tan/3 = 10 fi > 0. 
The solid (blue) lines mark the bounds on the cross sections allowed by 
the relic-density constraint 0.1 < tt x h 2 < 0.3 alone. For any given value of 
only a restricted range of mo is allowed. Therefore, only a limited 
range of mo, and hence only a limited range for the cross section, is allowed 
for any given value of m x . The thicknesses of the allowed regions are due 
in part to the assumed un certainti es in the nuclear inputs. These have 
been discussed at length in | 143 | 14 2| Q n ^ ne Q^gj. h anc j j a broad range of 
m x is allowed, when one takes into account the coannihilation 'tail' region 
at each tan/3 and the rapid-annihilation 'funnel' regions for tan/3 = 35, 50. 
The dashed (black) line displays the range allowed by the b — > sj constraint 
alone. In this broader range of mo and hence the spin-independent 

cross section is possible for any given value of m x . The impact of the con- 
straint due to rrih is shown by the dot-dashed (green) line. Comparing with 
the previous constraints, we see that a region at low m x is excluded by m/,, 
strengthening significantly the previous upper limit on the spin- independent 
cross section. Finally, the dotted (red) lines in Fig. 1 181 show the impact of 
the — 2 constraint. This imposes an upper bound on mi/2 and hence m x , 
and correspondingly a lower limit on the spin-independent cross section. 
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b igure 19. Limits from the CDM^HHland EdelweissESl experiments on the neutralino- 
proton elastic scattering cross section as a function of the neutralino mass. The Edelweiss 
limit is st ronge r at higher m x . These results nearly exclude the shaded region observed 
by DAMv414Zl The theoretical predictions lie at lower values of the cross section. 

This analysis is extended in pa nel ( b) of Fig. ^| to all the values 8 < 
tan (3 < 55 and we find overall that I I 

2 x 1(T 10 pb < a S i < 6 x 1CT 8 pb, (60) 
2 x 1(T 7 pb < a SD £ 10" 5 pb, (61) 

for fi > 0. (<tsd is the spin-dependent cross-section not shown in the 
figures presented here.) As we see in panel (b) of Fig. 1181 provides 
the most important upper limit on the cross sections for tan/? < 23, and 
b — ► S7 for larger tan /3, with — 2 always providing a more stringent 
lower limit than the relic-density constraint. The relic density constraint 
shown is evaluated at the endpoint of the coannihilation region. At large 
tan/?, the Higgs funnels or the focus-point regions have not been considered, 
as the ir lo cations are very sensitive to input parameters and calculational 
detaiira 

The results from a CMSSM and MSSM analysiJ23HII for tan (3 = 3 
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and 10 are compared with the most recent CDMePSI and Edelweis^^SI 
bounds in Fig. 1191 These results have nearly entirely excluded the regio n 
purported by the DAMPISH experiment. The CMS SM p redictioifSSl j s 
shown by the dark shaded region, while the NUHM cas o lis shown by the 
larger lighter shaded region. Other CMSSM results^^ are also available. 

I conclude by showing th e pro spects for direct detection for the bench- 
mark points discussed aboveSSI. Fig. [2U| shows rates for the elastic spin- 
independent scattering o f sup ersymmetric relics , including the proje cted 
sensitivities for CDMS iffSOl aim CRESSlEU (solid) and GENIUS^ 
(dashed). Also shown are the cross sections calculated in the proposed 
benchmark scenarios d iscus sed in the previous section, which are consider- 
ably below the DAM^ffSI 

range (10 5 — 10 6 pb). Indirect searches for 
supersymmetric dark matter via the products of annihilations in the galac- 
tic halo or in side the Sun also have prospects in some of the benchmark 
scenario 
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Figure 20. Elastic spin-indep ende nt scattering of supersymmetric relics on protons cal- 
culated in benchma rk sce nario compared with the projected sensitivities for CDMS 
IlE01 arid CRESSlEII (soli d! an d GENIU^Sl (dashed) . The predictions of our code 
(blue crosses) and Neutdrivei (red circles) for neutralino-nucleon scattering are com- 
pared. The labels A, B, ...,L correspond to the benchmark points as shown in Fig. 1121 
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